Improved Nearside-Farside Decomposition of Elastic Scattering Amplitudes 
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scattering matrix element Si by: 

ai = {2l + l){Si-l). 



We also recall that the PWS in (|l]), considered as a dis- 
tribution, is convergent if Si is asymptotically Coulombic 

i- 

The Fuller NF decomposition is realized by splitting 
P/(cos6'), considered as a standing angular wave, into 
traveling angular wave components 



A simple technique is described, that provides improved nearside-farside (NF) decompositions of 
elastic scattering amplitudes. The technique, involving the resummation of a Legendre partial wave 
series, reduces the importance of unphysical contributions to NF subamplitudes, which can arise in 
more conventional NF decompositions. Applications are made to a strong absorption model and to 
^ 16q _,_ i2(-, op^jgg^i potential at Siab = 132 MeV. 

PACS numbers: 24.10.Ht, 25.70.Bc, 34.50.-s, 03.65. Sq 

In heavy-ion, atomic and molecular collisions, an elas- 
tic differential cross section cr(0), where 6 is the scattering 
angle, is often characterized by a complicated interference 
pattern. This complicated structure makes it difficult to 
understand the physical phenomena involved in the scat- 
tering process, as well as the links between a{0) and the 
properties of the model that describes the phenomenon. 

In some cases, semiclassical methods explain the 
scattering pattern as the interference between simpler, 
and slowly varying, subamplitudes. If we ignore the com- 
plication that in some angular regions uniform asymp- 
totic techniques are necessary, then the semiclassical sub- 
amplitudes arise mathematically from saddle points or 
poles which account physically for the contributions from 
reflected, refracted or generalized diffracted semiclassical 
trajectories Q. The subamplitudes can be conveniently 
grouped into two types: those arising from semiclassical 
trajectories which initially move in the same half plane as 
the detector (N or nearside trajectories) and those from 
the opposite half plane (F or farside trajectories). 

Semiclassical methods are not always simple to apply 
and sometimes they have a limited range of applicability. 
Their limitations are determined by the range of validity 
for (presently known) asymptotic techniques that pro- 
vide an accurate approximation to the original quantum 
mechanical problem. 

In order to overcome these difficulties. Fuller]^ pro- 
posed more than 25 years ago, a simple NF decomposi- 
tion of the elastic scattering amplitude f{9). We write 
f{9) in the usual way as a partial wave series (PWS) 



(2) 



where (for 7^ 0, tt) 



P/ (cos 61) = q\ ^ (cos 61) + q|+^ (cos 61), (3) 
Qp^(cos^) = ]-[Pi{cose) ± -Q,(cos^)], (4) 

Z TT 

with Qi{cos9) the Legendre function of the second kind 
of degree I. 

Inserting (||) into ([^), splits f{9) into the sum of two 

subamplitudes f'-^HO). For ^ sin 61 > 1, the (cos 61) 
behave as 



q|^^ (cos 61) 



1 



2nX sin 9 



exp[Ti(A0--)], (5) 



with X = I 



This asymptotic behavior suggests that 
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where k is the wavenumber, Pi{cos9) is the Legendre 
polynomial of degree I and a; is given in terms of the 
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the f^^\9) should correspond to NF trajectories respec- 
tively appearing in the complete semiclassical decompo- 
sition of f{9). 

However the NF decomposition is in general less sat- 
isfactory than the full semiclassical one. For example, if 
two or more N or F semiclassical trajectories contribute 
to the same 9, interference effects may appear in the N or 
F cross sections. The Fuller NF decomposition has, how- 
ever, the merit of being simple and, although inspired by 
the semiclassical theories, it uses only quantities calcu- 
lated within the exact quantum mechanical treatment. 
The NF method therefore bypasses problems associated 
with the applicability and validity of semiclassical tech- 
niques. 

The physical meaning attributed to the f^^\9) is im- 
plicitly based on the (unproven) assumption that it is 
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possible to perform on the PWS, written in terms of the 
q[^\cos0), the same manipulations which are used in 
the complete semiclassical decomposition of f{0). These 
manipulations are essentially path deformations in A of 
the integrals into which (|l|) can be transformed, using ei- 
ther the Poisson summation formula or the Watson trans- 
formation. The consequences of these path deformations 
depend on the properties of the terms in the PWS when 
they are continued to real or complex values of A from 
the initial half integer A values. The splitting of Pi{cos9) 

into q[^\cosO) modifies these properties and can cause 
the appearance of unphysical contributions in the f^^\0) 
which cancel out in f{9). 

In spite of these possible limitations, the Fuller NF 
decomposition is widely used. In many cases ||, it de- 
composes f{9) into simpler subamplitudes which are free 
from the unphysical contributions that can arise from 
the above mathematical difficulties. However for a few 
examples, the NF subamplitudes can be directly com- 
pared with the corresponding semiclassical results and 
it is found that the Fuller and semiclassical decompo- 
sitions predict different results. One classic example is 
pure Coulomb scattering. For repulsive Coulomb poten- 
tials only a N contribution is expected semiclassically ( jl] , 
p. 56), whereas the Fuller NF decomposition yields also 
a F contribution . Another example is the angular dis- 
tribution for a strong absorption model (SAM) with a 
two parameter (A and A) symmetric S matrix element 
and Fermi- like form factors Q . For a fixed value of the 
cut-off parameter A and for a sufficiently large value of 
the diffuseness parameter A, the Fuller NF cross sections 
agree with the semiclassical results only up to a certain 
value of 9, which decreases with increasing A. 

Fortunately, the Fuller NF subamplitudes contain in- 
formation which allows one to recognize the unphysical 
nature of the undcsired contributions. Suppose f^^\9), 
or H^)) contains a single contribution from a station- 
ary phase point at X{9). Then the derivative with respect 
to 9 of the phase of /(+'(6'), or f^-^9\ is equal to \{9), 
or -X{9) respectively; it depends on 9 p. 57). Follow- 
ing Fuller we will call this derivative the Local Angular 
Momentum (LAM) for the N (or F) subamplitude. Only 
for certain generalized diffracted trajectories is the LAM 
expected to be constant, equal to the angular momentum 
of the incoming particle responsible for the diffraction. In 
the semiclassical regime, this constant value is expected 
to be large. Because of this, if we observe in a certain 9 
range that LAM « 0, this can be considered the signature 
of the unphysical nature of the N or F subamplitudes in 
that range of 9. This occurs for the LAM of the Fuller 
Coulomb F subamplitude, and for the NF subamplitudes 
of the SAM in the angular region where the NF cross sec- 
tions differ from the semiclassical results. In both cases 
this decoupling of 9 from LAM suggests the unphysical 
nature of the subamplitudes. Thus an analysis of the 
LAM can avoid misleading interpretations of cross sec- 
tions obtained from the Fuller NF decomposition. 

However the problem of obtaining more satisfactory 



NF decompositions remains open. A possible solution 
to the problem was proposed by Hatchell 0| , who used a 
modified NF decomposition. The modifications consisted 
of, first, in writing f{9) in the resummed form {9 ^ 0) 
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P/(cos^ 



(6) 



r — 1,2,..., and, second, in using a different splitting 
for the Legendre polynomials into traveling waves. 

The use of the resummed form (|6|) for was orig- 
inally proposed ^ by Yennie, Ravenhall, and Wilson 
(YRW) to speed up the convergence of the PWS for high- 
energy electron-nucleus scattering. Equation ^ is an 
exact resummation formula, of order r, which is derived 
from the recurrence relation for Legendre polynomials. 
The YRW resummation formula can be derived by iter- 



ating r times, starting from 
identity 
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a;, the resummation 



^a«PKcos^?), (7) 

i=0 



where 
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with a_i ~ 0. 

Note that does not depend on r, unlike the 

Fuller NF subamplitudes which do depend on the value 
of r used. This is a consequence of the property 
lQi-\{co%9) — > 1 as Z ^ jo). In the Hatchell approach, 
the dependence on r arises because the functions used in 
place of the Qi(cos0) obey a recurrence (inhomogencous) 
relation different from that for (cos 61). It is worth- 
while to note that (^ for r > 1 allows one to drop, for 
6* > 0, the 1 appearing in the term Si ~ \ in (|2|). Fur- 
thermore for r > 1, (^ produces a convergent PWS even 
when Si is asymptotically Coulombic JToj . 

Using his method, Hatchell has shown j?) that the un- 
physical contributions to the SAM NF cross sections sys - 
tematically decrease on increasing r. More recently 

it was shown that, using the Fuller qI^"* (cos 6') func- 
tions in (PI), gives even better results. The superiority 



of the (5|^''(cos0) functions seems to be connected with 

the greater rapidity with which the oj^' (cos 0) approach 
their asymptotic behavior (|^) compared to the 

Hatchell NF functions. 

The success of using (||) before applying the NF decom- 
position depends upon the properties of the a) . For 



the SAM, the contributions from low / values rapidly 
decrease with increasing r. As a result, the most im- 
portant partial waves move to higher values of I, where 
a semiclassical description is physically more reasonable. 

However, in some cases, (||) acts in the opposite di- 
rection, by enhancing the undesired unphysical contribu- 
tions to the NF subamplitudes. We have found that this 
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happens, for example, for pure Coulomb scattering, for 
scattering by an impenetrable sphere, and for the SAM 
(see ) when the cross section is calculated at an angle 
TT — 9, using the property P; [cos(7r — 6*)] = (— l)'P/(cos 6*). 

Onepossible solution to this intriguing puzzle is to re- 
gard (0) as a particular case of the modified resummation 
identity @ 

oo ^ oo 

J2'^t'^Pi{cose) = -^a«P,(cos0), (9) 

'^-^ ai + Pi cos H ^-^ 



1=0 



1=0 



with ai + Pi cos 9 and 



For a;, j3i ^ 0, the r.h.s. of depends only on the ratio 
Pi/ai. Thus, without loss of generality, we can assume 
ai = 1 for all i. By iterating (j^) r times, we can write 
f{0) in the modified resummed form 



I 



(i-l) , 



Z + 1 




^ar^Pi (cos 0), (11) 



r = 1,2,.... The YRW resummation formula (||) is ob- 
tained with P\ — Pi = . . . = P,. — —\. 

The resummation identity (^) is a particular case of a 
more general one |l^, which uses a basis set of reduced 
rotation matrix elements; this gives the amplitude for 
more general scattering processes than those described by 
(|l|) . For these general PWS, a Fuller- like NF decomposi- 
tion can be introduced |l^, |l^, ^ which allows the scat- 
tering amplitude to be split into NF subamplitudes. In 
some cases, the NF cross sections contained unexpected 
(unphysical) oscillations |]l^ , which are enhanced if the 
generalization of (^) is used, but which disappear for an 
appropriate choice of the /3-parameter in the generaliza- 
tion of (^. 

The considerable successes achieved by the original 
Fuller NF decomposition suggests that the modified re- 
summed form ( |ll|) be used to diminish unphysical con- 
taminations to the NF subamplitudes when they are 
present. To do this, we must give a practical rule to 
fix the values of the /3-parameters. In Refs. ^ it 
was proposed to select the value of /3 = /3i = . . . = /?r 
so that (1 -I- /3cos0)~'' approximately mimics the shape 
of the angular distribution. The shape of the cross sec- 
tion can however be very different from that given by 

+ P cos 9)~^. It is therefore desirable to test a different 
recipe, based on a simple rule. The quantitative recipe 
proposed here is inspired by the observation that the 
modified resummation formulas produce a more physi- 
cal NF decomposition by reducing the contribution from 
the low I values in the resummed PWS. This suggests 
we select the /32, • • ■ , /3r in r repeated applications of 
(^ , so as to eliminate as many low partial waves as pos- 
sible from the final PWS in (|Tl|). The transformation 



from {a[* to is linear tridiagonal, with coeffi- 

cients linear in Pi, which means application of r succes- 
sive resummations allows one to equate to zero the lead- 

(r) 

ing r coefficients al , with / = 0, 1, . . . , r — 1, by solving 
a system of r equations of degree r in the parameters 
Pi,P2, ■ ■ ■ ,Pr- We will call the resummation defined in 
this way an improved resummation of order r. 

It is straightforward to show that the improved resum- 
mation of order 1 is obtained by choosing 

Pi = -3ao/ai, (12) 

while the improved resummation of order 2 is given by 

Pi,2 - (P ± VB^ - 4A)/2 , (13) 
with A and B solutions of the linear equations 



{lao + ^a2)A + 



-a-o 



{^ai + ^a3)A + (ao + fa2)P = -ai. 



(14) 



Higher order improved resummations require the solution 
of more complicated systems of equations. 

In all the cases we have analyzed using r < 2, we find 
that the improved resummations considerably reduce the 
width of the angular regions in which the Fuller NF cross 
sections exhibit unphysical behavior. In these analyses, 
we have used Si from simple parametrizations as well as 
from some of the optical potentials currently employed to 
describe light heavy-ion scattering. We show below our 
results for two particular examples. The first example is 
a SAM (Fig. 1), whilst the second example is the ^^O -I- 
i^C coUision, at Piab = 132 MeV, using the WSl optical 
potential of Rcf. p9|] (Fig. 2). For both these cases we 
have dropped the 1 in the term S'; — 1 in Eq. (||). The 
calculations were performed applying: first, an improved 
resummation of order r = 0, 1, 2, with r = meaning no 
resummation, second, the Fuller NF decomposition (|^), 
(^) , and third a YRW resummation of the NF subampli- 
tudes using the extension of (|^) to the linear combination 
(^) of integer degree Legcndrc functions of the first and 
second kinds [^. This latter resummation ensures the 
convergence of the final NF PWS. The results obtained 
from these three steps will be indicated by the notation 
R = Oy, 1y, 2y. 

For the SAM we have chosen the parameters to be 
A = 10 and A = 2. In Fig. 1 we have plotted the di- 
mensionless quantity Ak'^a{9) sinO since the correspond- 
ing NF semiclassical quantities are expected to have a 
pure exponential slope |Q. Furthermore, because the Si 
are real, f{9) has a constant phase (and its phase deriva- 
tive is of no interest), and the /'•^H^) have the same 
moduli but opposite phases. Thus we need only show 
the N, or F, LAM and similarly for the cross sections. 

Figure 1 also shows the results obtained on introduc- 
ing the NF decomposition directly into (|l|) and without 
dropping the 1 in (||). The NF subamplitudes obtained 
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FIG. 1: Strong absorption model N (continuous lines) and F 
(dashed lines) cross sections (lower panels) and LAM (upper 
panels) calculated using the R = 0,0y,1y,2y NF decompo- 
sitions. The thin curve shows the full cross section. The 
thin dotted curve shows the F {— N) cross section (displaced 
downward by one unit) for the unphysical amplitude /j^' (6) . 



in this way are rapidly convergent and no final YRW re- 
sumniation is needed. These results are indicated by the 
notation R = 0. 

For the original Fuller NF method, R = 0, and for 
the case R = ly, we have plotted the F cross sections 
and LAMs (dashed curves); for the cases R = Oy and 
R = 2y the N quantities (continuous curves) are dis- 
played. Figure 1 shows that for R = the unphysical 
contributions dominate the F (= N) cross section over 
most of the angular range. The expected exponential be- 
havior is not present in the F cross section curve and the 
F (= -N) LAM « for 6* > 60°. However, at smaller 
angles, oscillations in the F LAM curve indicate that an- 
other contribution is present which interferes with the 
unphysical one. This behavior does not support the con- 
jecture that the F LAM of this other contribution has 
the semiclassical value A. The major part of the unphys- 
ical contribution, which dominates the F subamplitude, 
is fs'^\d) = -[27rfc(l - cose)]-\ This is the F compo- 
nent of fs{0) = iS{l — cos6)/k, obtained by dropping Si 
in the term Si — 1 in (||). The cross section of fg^\o) 
is shown, downward shifted by one vertical unit, by the 
thin dotted curve in Fig. 1. 

The R = Oy method provides more satisfactory results, 
which are rather good at forward angles. With the ex- 
clusion of a small region around = 0°, where (||) does 
not hold, and up to 6* ss 50°, the N (= -F) LAM agrees 
closely with the expected semiclassical value of —A and 
the N (= F) cross section curve follows the expected ex- 
ponential behavior. For 9 > 120°, the N cross section is 
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FIG. 2: Optical potential model N (continuous lines) and F 
(dashed lines) cross sections (lower panels) and LIP (upper 
panels) calculated using the R = Oy and R = ly NF decom- 
positions. The thin curves show the cross section and LIP 
obtained using the full quantum amplitude in the left panels, 
and the N (continuous) and F (dashed) cross sections and LIP 
using classical mechanics in the right panels. The indices in 
the right panels identify the curves corresponding to different 
branches of the classical deflection function. 



still dominated by an unphysical contribution. At inter- 
mediate angles, 50° ^ < 120°, interference oscillations 
appear both in the N cross section and in the N LAM 
curve. It is interesting to note that the LAM is more 
sensitive to interference effects than is the cross section. 
Also, in the interference region, one cannot attach the 
meaning of a local angular momentum to the subampli- 
tude phase derivative. In our case, in this interference 
region, the N LAM curve oscillates around the expected 
semiclassical value of —A in the region, 50° d < 80°, 
where the true semiclassical component dominates the N 
subamplitude, and around the unphysical value of at 
other angles. 

The effectiveness of the improved resummation proce- 
dure is evident in the right panels of Fig. 1. Using the 
R = ly method ( for which Pi = -0.800) the F (= - 
N) LAM and the F (= N) cross section are in agreement 
with the semiclassical results up to 6* « 120°. For R = 2y 
(which has /3i.2 = —0.879 ± 0.076 i) the agreement covers 
almost the whole angular range. The small irregular os- 
cillations appearing at large 6 for the N LAM curve, with 
R = 2y, probably arise from the precision limitations (64 
bit floating point representation) of the calculations. 

Figure 2 shows our results for the optical potential. In 
the upper panels we display LAM/fc, which we call the 
Local Impact Parameter (LIP), and in the lower panels 
a Log plot of a{6) sin 6. The left panels show the results 
for the usual Fuller NF decomposition, R = Oy. The thin 
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continuous lines, in the left panels, show the cross section 
and LIP for the full amplitude. 

The behavior of the NF LIP curves is mostly simpler 
than that of the fuU amplitude LIP. At 6* w 90°, the 
R = Oy N LIP curve oscillates around 0, indicating the 
possible dominance of an unphysical contribution. This 
contribution also appears to be responsible for oscilla- 
tions in the N LIP curve around other values (different 
from 0), and for oscillations in the N cross section for 
9 > 30°. These oscillations are absent in the N curves 
in the right panels where the results for R = ly are 
shown (/3i = -0.999 - 0.099 i). Both the N cross section 
and N LIP curves for R = ly are considerably simpler 
than those obtained using R = Oy, while the F curves 
are essentially the same; an exception is the less oscil- 
latory F LIP for > 120°. This indicates that, apart 
from 9 > 120°, the unphysical contribution for R = Oy 
has a modulus which is much smaller then that of the 
F semiclassical subamplitude. We have also applied the 
improved resummation R = 2y. The results are practi- 
cally the same as those obtained using R = ly and are 
not shown. 

The cleaning by the R = ly procedure of the origi- 
nal R = Oy NF subamplitudes is impressive and allows a 
clear identification, in the NF cross sections at 9 > 120°, 
of the dominance of semiclassical trajectories refracted 
from the nuclear part of the interaction. In the right 
panels of Fig. 2, this interpretation is confirmed by the 
agreement, for 9 > 120°, between the NF curves and the 
corresponding classical mechanical results (thin lines 1 
and 2). The thin lines show, in the upper panel, differ- 
ent NF branches of the the impact parameter and their 
dependence on 9 (with appropriate signs) using only the 



real part of the complete interaction. In the lower panel 
we show the classical contributions to the cross section 
from these branches, in which we have included in the 
usual simple way ([|l[, p. 49) the absorptive effects of the 
imaginary part of the optical potential. 

Our new resummation NF procedure clearly improves 
the original Fuller NF decomposition, as is evident in 
the examples presented here. On the one hand, our re- 
sults confirm the importance of NF decompositions for 
gaining insight into the properties of the subamplitudes 
responsible for complicated structures in cross sections. 
On the other hand, they confirm the empirical origin of 
NF decompositions and suggest caution in the interpre- 
tation of results obtained from NF techniques. However, 
different NF decompositions can be used to check what 
parts of the resulting NF subamplitudes are independent 
of the particular technique used. Only properties stable 
with respect to different NF decompositions, can be con- 
sidered as manifestation of some physical phenomenon. 
In addition, we have shown that it is desirable to inves- 
tigate the behavior of the LAM. This quantity is more 
sensitive to interference effects than are the NF cross sec- 
tions, and a null value (or an oscillatory behavior around 
zero) of the LAM in a certain angular range indicates an 
unphysical contribution. 
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